Random and aperiodic quantum spin chains: A comparative study 
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According to the Harris-Luck criterion the relevance of a fluctuating interaction at the critical point 
is connected to the value of the fluctuation exponent u). Here we consider different types of relevant 
fluctuations in the quantum fsing chain and investigate the universality class of the models. At the 
critical point the random and aperiodic systems behave similarly, due to the same type of extreme 
broad distribution of the energy scales at low energies. The critical exponents of some averaged 
quantities are found to be a universal function of u), but some others do depend on other parameters 
of the distribution of the couplings. In the off-critical region there is an important difference between 
the two systems: there are no Griffiths singularities in aperiodic models. 
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The presence of quenched disorder has a profound ef- 
fect on quantum phase transitions, especially in onc- 
dimensionsEJ. A common feature of these one-dimensional 
random magnets is the extremly broad distribution of en- 
ergy scales at low energies in the vicinity of the critical 
point. This property makes it possible to use approxi- 
mate renormalization-group transformationsB , which be- 
comes exact at the critical point. Recently, Fishertl have 
obtained new striking results for the random transverse- 
field Ising spin chain defined by the Hamiltonian: 
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Here the a'f , of are Pauli matrices at site I and the Ji ex- 
change couplings and the hi transverse-fields are random 
variables. The random couplings (and fields) in (Q) show 
unbounded fluctuations, the cummulated deviation from 
the average coupling [J] av grows asymptotically with the 
size L as: 



A(L) = ( Jl ~ Wav) ~ L" 



(2) 
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with the fluctuating or wandering exponent u> — cj raru j = 
1/2. A similar type of unbounded fluctuations occur in 
aperiodic sequences, too, which are generated through 
substitutional rules. For example the Rudin-Shapiro 
(RS) sequence is built on four letters A,B,C and D with 
the substitutional rule: 



AB , B^AC 



DB 



D^DC, (3) 



thus starting with a letter A one proceeds as: A — > 
AB -> ABAC -> ABACABDB -> etc, and one may 
assign different couplings to the different letters. The 
wandering exponent of the sequence 



ln|A 2 |/lnA 1 



(4) 



which is given in terms of the leading and the next-to- 
lcading eigenvalues Ai : 2 of the substitution matrixfl, is 
^RS = 1/2, i.e. just the same as for the random sequence. 



Then naturally the question arises, which type of crit- 
ical behavior can be found in quantum Ising chains with 
an unbounded aperiodic modulation in the couplings (or 
fields). It is knowjjp.from the Harris-Luck relevance- 
irrelevance criterionEnj that the phase transition in an 
aperiodic (or quasi-periodic) quantum Ising chain be- 
longs to the Onsager universality class, only if the fluctu- 
ations in the couplings are bounded, uj < 0. For marginal 
sequences with uj = non-universal critical behavior and 
coupling dependent anisotropy exponent was found in ex- 
act calculationsEra. 

In the relevant situation with uj > there are only 
a few exact result^, however, on the basis of scaling 
considerationsBliiJiij'Ej the phase transition in relevantly 
aperiodic and random systems seem to be similar. For 
instance the energy gap of a finite, aperiodic chain at the 
critical point scales as: 



AE(L) ~ exp(-const • W) 



(5) 

just as in the random modelHOIi3. Thus one could ask 
the question, whether the wandering exponent ui alone is 
sufficient to characterize the universality classes of mod- 
els with fluctuating interactions; and whether or not the 
critical exponents of the random and the RS chains are 
the same? 

To answer to these questions in this Letter we are going 
to study systematically the critical behavior of aperiodic 
quantum spin chains with unbounded fluctuations and 
compare it with that of the random system. We mainly 
concentrate on the RS model, but some results are also 
presented for a family of relevant sequences defined on 
k- letters (fc-general sequence) Ai, Ai,..., Ak, with the 
substitutional rules Ai — > Ai_iAj + i for i < k/2 and 
Ai — > Aj + iAi_i for i > k/2 (one identifies Ao = Ai 
and Ak+i = Ak). The wandering exponent (f|) for the 
fc-general sequences, given by 



Lu k = m[2cos(7r/fc)]/ln2 , 



(6) 



is positive for k > 4. It is easy to see that the fc-general 
sequence for k — 2 and k — 4 is just the Thue-Morse and 
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the Rudin-Shapiro sequence, respectively, and obviously 

^4 = ^RS = 1/2- 

In the Hamiltonian ([!]) we have chosen homogeneous 
fields hi = h and two- valued couplings: Ji — X for letters 
Ai with i < (k + l)/2, and J; = 1/A for letters with i > 
(k + l)/2. For an odd k we take J; = 1 for i = (k + l)/2. 
Similarly, for the random model we take homogeneous 
fields and binary distribution of the couplings: J; = A 
and J i — 1/A, with th^^ame probability. The critical 
point, simply given by 



[In J] av = [In h] 8 



(7) 



is at ho = 1, both for the random and aperiodic systems. 

In the following we calculate different physical quan- 
tities (energy-gap, surface and bulk magnetization, crit- 
ical magnetization profiles, correlation length, etc.) for 
the RS-chain and compare with the known results on the 
random chain. We also present some results for the k- 
general chains, more details will be given in a subsequent 
publicationE^I. In the actual calculation, we first trans- 
form H in (Q) into a fermionic modelt3 and then use the 
representation iro, which necessities only the diagonal- 
ization of a 2L x 2L tridiagonal matrix. The energy gap 
then corresponds to the smallest positive eigenvalue of 
the matrix, whereas the magnetization profiles can be 
obtained from the corresponding eigenvectors. Details 
will be presented elsewhereEj. 

We are going to calculate the average quantities of fi- 
nite systems both for the random and aperiodic chains. 
For the aperiodic systems the average is performed as 
follows. We generate an infinite sequence through sub- 
stitutions, cut out sequels of length L, which start at 
all different points of the sequence and then average over 
the realizations. Among these realizations there are R(L) 
different, which grows linearly with L: R(L) = ah. For 
the RS-chain the number of different realizations is less 
than 16L, and one can obtain the exact average value 
by the following procedure. Consider the RS-sequence in 
(||) generated from a letter A, take the first 4L sequels of 
length L, starting at positions 1, 2, . . . , 4L and take also 
their reflexion symmetric counterparts. Then repeat the 
procedure starting with a letter D. The averaging then 
should be performed over these 16L realizations. The 
fact that the number of different samples grows linearly 
makes it possible to obtain numerically exact average re- 
sults for relatively large (L < 512) aperiodic chains. On 
the other hand, for random chains with the binary cou- 
pling distribution there are 2 L different realizations and 
we performed the average over some (N = 50000) ran- 
domly chosen samples. 

First, we investigate the probability distribution of the 
energy gap Pl{AE) at the critical point of the RS-chain. 
According to eq([|) the appropriate scaling variable is 
InAE/L 1 ^ 2 . Indeed, as seen in Fig.l, the accumulated 
probability distribution function 
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1. Scaling plot of the accumulated probability density 



fi(ln AE) versus the scaling variable (In AE/L 1/2 ) for the RS 
sequence (exact average) with A = 4, c.f. eq.(^). The insert 
shows the same for the the random sequence (a; and y axis as 
in the main figure). 

with P l QxlAE) = P L (AE)AE, has a good data col- 
lapse using this reduced variable. Considering now the 
same quantitity for the random chain one can observe the 
same scaling behavior, as shown in the insert of Fig.l. 
Thus we can conclude, that both the random and the 
RS-chains have logarithmically broad distribution of the 
energy gaps at the critical point, from this fact one ex- 
pects similar consequences for the critical behavior in the 
two systems. 

Next we turn to consider the surface magnetization of 
the quantum. Ising chains, which can be obtained from 
the formulaO 



m s (L, h) = 



L-l I 

i + EIKt 

1=1 j=l 



-1/2 



(9) 



for a given realization. For one single RS-chain, start- 
ing with-,the letter A this expression has been evaluated 
exactlyiMi with the following results. For A > 1 there is a 
finite critical point surface magnetization, which goes to 
unity as A — > oo, whereas for A < 1 the surface magneti- 
zation vanishes as 



m s (L, ho) ~ exp(— const • ^/L) . 



(10) 



Now considering all the possible realizations of the chain 
at the aperiodic critical point eq(0) we will see, that 
in typical samples i.e. which occour with probability 
one, the surface magnetization vanishes as in cq(fL0|), 
whereas there are rare events with a vanishing proba- 
bility Prarc(i) ~ L _7 , but with a surface magnetization 
of order unity. Thus the critical surface magnetization 
is non- self averaging, it is determined by the rare events, 
and its scaling dimension x s m defined by the asymptotic 
relation [m s (L, h )] av ~ L~ x ™ is just x s m = 7. 

In the following, using this observation, we calculate 
x s m exactly. First, we note that to each sample with a 
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given distribution of the couplings one can assign a walk, 
which starts at the origin and takes the l-th step +1 
(— 1) for a coupling J; = A (J; = 1/A). Then, as shown 
in Reia, taking the limit A — * oo only those samples 
have non- vanishing surface magnetization, where the cor- 
responding walk never visits sites with negative coordi- 
nates. Thus the proportion of rare events is given by the 
surviving probability of the walk: P rare (L) = P SUIV (L). 
For the random chain one should consider the random 
walk with P surv llA^u L -1 / 2 , thus one immediately gets 
the exact resultBl-H 

< n (rand) = 1/2 . (11) 

For the RS-chain one can perform the exact analysisEl, 
from which here we present the leading finite-size de- 

pendence: [m s (L, h )} av = § + i) L' 1 / 2 + 0{L^/ i ) 

thus the corresponding scaling dimension 



,(RS) = 1/2 



(12) 



is the same as for the random model. 

For general relevant aperiodic sequences the surface 
magnetization scaling dimension can be obtained by the 
following scaling consideration about the surviving prob- 
ability of the corresponding aperiodic walk. Let us per- 
form a discrete scale transformation, which corresponds 
to a substitutional step of the sequence when the length 
of the walks scales as L — ► LAi, whereas the trans- 
verse fluctuations as y — > y|A 2 |. Then N(L) the num- 
ber of i-steps surviving walks from the total of different 
R(L) = aL walks scales as N(L) -> iV(Aii) = \A 2 \N(L), 
since the number of these walks is proportional to the size 
of transverse fluctuations. Consequently the surviving 
probability P BUIV (L) = N(L)/R(L) satisfies the scaling 
relation P surv (AiL) = |A2|/AiP surv (£), with the general 
result 



= 1-0) 



(13) 



This simple exponent relation, which is valid for the ran- 
dom and RS chains, has been numerically checked for the 
first members of the /c-general chain. 

Next we turn to calculate the correlation length critical 
exponent v from the 6 = In h dependence of the surface 
magnetization. In the scaling limit L ^> 1, \6\ <C 1 the 
surface magnetization can be written as [m s (L, 5)] av — 
[m s (L,0)] m m s (5L 1 '' / ). Expanding the scaling function 
into a Taylor series m s (z) = 1 + bz + 0(z 2 ) one obtains 
for the (5 correction to the surface magnetization: 



[m s (L,6)] av 

[m s (L, 0)] 



(14) 



with = Xjv — x s m . This exponent can also be deter- 
mined ex^cidy in the A — > oo._Limit from random walk 



argumentaifl. As shown in Reaia the surface magnetiza- 



tion of rare events is given by: 
m s (L,S) = (1 + n)- 1 / 2 -5 



EI 



(n + l) 3 /2 



+ 0(6 2 ), (15) 



where the corresponding surviving walk returns n-times 
to its starting point after li, 1%, ■ ■ ■ , l n steps. On the basis 
of (|l^) we argue that is connected to the surface fractal 
dimension d s of the surviving walks, defined through the 
asymptotic dependence of the number of return points n 
of a surviving walk of L steps, n ~ L ds . Now, to perform 
the average of the linear term in ([l5]) one should take into 
account that from the R{L) different samples there are 
0(n), which give the most important contribution: each 
of those has 0(n) return points of characteristic lengths 
li ~ L. Consequently, the average of the linear term in 
( |l5| ) grows as L ds ^ 2 , thus comparing with (14) one gets 
the exponent relation: 
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(16) 



The random case ^(rand) = 2 is formally contained in 
( |l6| ) with d s = 0, since a surviving random walk returns 
n = 0(l)-times to the starting point. For the family of 
fc-general sequences d s = 1/2 for all values of k, thus 
the corresponding correlation length exponent is v{k) = 
4/(5 — 4cjfe), with Uj~ given in (^|). In particular for the 
RS-sequence we get 



v{RS) = 4/3 , 



(17) 



which we also checked numerically by evaluating eq. jl5| ) 
exactly up to L = 2 21 . Thus we conclude that the RS 
and the random Ising quantum chain are not in the same 
universality class. 

Next we turn to study the magnetization in the bulk 
and calculate the scaling dimension of the bulk magneti- 
zation x m — (3 /v from the behavior of the average mag- 
netization profiles at the critical point in fkiite systems. 
As will be shown in the detailed publicationtj the average 
magnetization profiles for the RS-chap-jajre in excellent 
agreement with the conformal result a^O. in a similar 
way as observed for random chainscSEia From the scal- 
ing behavior of the profiles be estimated as 



( (RS) = 0.160 ± .005 



(18) 



which also differs from the random chain value 
x m (rand) = (3 - \/5)/4 = .1910. 

Finally, we turn to discuss the properties of the RS- 
chain in the disordered phase. For the random model, as 
known from exact resultsBta, the susceptibility (diverges 
in a whole region, in the so called Grifnths-phasecJ. This 
non-analytical behavior is due to a broad distribution 
of the energy gaps, which scales in finite systems 
AE(L, 5) ~ L~ z ( s \ with the dynamical exponent z(5)\ 
As a consequence, in the infinite chain there is no fi- 
nite time-scale and the autocorrelation function decajffi 
algabraically with an exponent which is related to z((5)E3. 

For relevantly aperiodic chains the same type of 
scenario, i,e. the existence of a Griffiths-phase is 
speculated!! To clarify this issue for the RS-chain we 
have calculated the distribution function of the energy 
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FIG. 2. The accumulated probability density Q(AE) for 
the RS sequence (exact average) with A = 4 slightly above 
the critical point (h = 1.5). The distribution is chopped off 
at lnAE~l n (h = 1.5) ~ 7.3. The insert shows In AE~l n (h) 
versus the distance from the critical point (ft — ho ~ 8 for 
h — > /io) in the RS sequence for different system sizes. The 
srtraight line has slope —1, as predicted by nWv. 



gap in the disordered phase. As seen on Fig. 2 the ac- 
cumulated distribution function for the RS-chain has 
a qualitatively different behavior, than for the random 
chain: there is an L- independent cut-off at AE m i n (5). As 
a consequence in the disordered phase the autocorrelation 
function always decays exponentially with a relaxation 
time t r ~ AE~l n (S) and the susceptibility and other 
physical quantities remain analytic. Therefore there is no 
Griffiths-region in the RS-chain and a similar scenario is 
expected to hold for any aperiodic quantum spin chain. 

One can estimate the minimal energy gap AE m i n (S) 
close to critical point from the related expression in eq(|^) 
by replacing there L with the corresponding characteris- 
tic length-scale in the off-critical region. We note, that 
this length-scale is not the average correlation length, 
since it is related to a single sample.. .In this case the 
aperiodic length scale Z ap is obtaincdOEj by equating the 
aperiodic (fluctuating) energy contribution ~ l^ 1 with 
the "thermal" energy ~ S giving Z ap ~ S^ l /^ x ~ w \ Thus 
the minimal energy gap is expected to scale as 



AE m i n (5) ~ exp(— const • S 



(19) 



This relation is indeed satisfied for the RS-chain, since 
in this case according to numerical results In A_E m i n (6) ~ 
<5 _1 , as can be seen in the insert of Fig. 2. 

To summarize, we have made a comparate study of 
the critical properties of random and relevantly aperiodic 
quantum Ising chains. At the critical point the two sys- 
tems behave very similarly, some critical exponents are 
found even identical. This similarity is mainly attributed 
to the fact that, at the critical point the distribution func- 
tions of different physical quantities in the two systems 
qualitatively agree. Outside the critical point, where the 
low energy tail of the distributions are of importance, the 



two systems behave differently. In the random system, 
with exponentially many realizations, there is no minimal 
energy scale, whereas in the aperiodic chain with linearly 
many independent realizations there exists a minimal en- 
ergy cut-off. This leads then to the absence of the Grif- 
fiths phase in aperiodic systems. 
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